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Abstract 

A stable neutralino x?) assumed to be the hghtest supersymmetric particle, is a 
favored particle physics candidate for the cosmological Dark Matter. We study co- 
annihilation of the lightest neutralino with the lighter scalar top quark ti. We show 
that for natural values of the neutralino mass, <300 GeV, the Xi"ii mass difference 
has to exceed ~ 10 to 30 GeV if Xi is to contribute significantly to the Dark Matter. 
Scenarios with smaller mass splitting, where ti is quite difficult to detect at collider 
experiments, are thus cosmologically disfavored. On the other hand, for small ti^Xi 
mass splitting, we show that co-annihilation allows very large neutralino masses, 

m^o ~ 5 TeV, without "overdosing" the Universe. 
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There is convincing evidence |T[] that most matter in the Universe is dark (non- 
luminous) , 

o.2<r]DM<i (1) 

where ^dm is the Dark Matter (DM) density in units of the critical density, so that Q = 1 
corresponds to a flat Universe. On the other hand, analyses of Big Bang nucleosynthesis 
[0] imply that most DM is non-baryonic (although dark baryons probably exist as well). 

One of the favorite particle physics candidates for DM is the lightest neutralino Xi 
[^, assumed to be the lightest supersymmetric particle (LSP). It is stable if R-parity 
is conserved this is also a sufficient (although not necessary) condition for avoiding 
very fast nucleon decay in supersymmetric theories. The LSP makes an attractive DM 
candidate since the primary motivation for its introduction comes from particle physics 
arguments supersymmetry stabilizes the huge hierarchy between the weak and Grand 
Unification scales against radiative corrections, and if it is broken at a sufficiently high 
scale, it allows to understand the origin of the hierarchy in terms of radiative breaking of 
the Standard Model (SM) electroweak SUl(2) x U(1)y gauge symmetry; furthermore it 
allows for a consistent unification of the gauge couplings. 

Supersymmetric contributions to DM then come as extra bonus, and for wide regions 
of parameter space, the LSP relic density falls in the preferred range eq.([l|). This is true 
in particular if the LSP is mostly a superpartner of the U(1)y gauge boson, i.e. bino- 
like, and if both m^o and the masses of SU(2) singlet scalar leptons fall in the natural 
range below a few hundred GeV |^ (but above ^ the mass range excluded by the LEP 
experiments). 

The previous statement assumes that XiXi annihilation reactions are the only processes 
that change the number of superparticles at temperatures around Tp ~ m^o/20, where 
the neutralino Xi decouples from the plasma of SM particles. It has been known for some 
time 1^ that this is not true if the mass splitting between the LSP and the next-to-lightest 
supersymmetric particle P is small. In this case, reactions of the type 

XS + X^P + F, (2) 

where X, Y are SM particles, occur much more frequently at a temperature T Tp than 
XiXi annihilation reactions do. The rate of the latter kind of process is proportional to two 
powers of the Boltzmann factor exp(— m^o/Tp) ~ exp(— 20), whereas for m^o ~ mp the 
rate for reaction (H) is linear in this factor. These reactions will therefore maintain relative 
equilibrium between the states Xi and P until long after all superparticles decouple from 
the Standard Model plasma. 

The total number of superparticles can then not only be changed by XiXi annihilation, 
but also by the "co-annihilation" processes 

X? + P ^ X + F and P + P^*^ ^ X + Y. (3) 

Eventually all particles P and P* will decay into Xi (plus SM particles). In order to 
compute today's LSP relic density, we therefore only have to solve the Boltzmann equation 
for the sum nsusY of densities Ui of all relevant species of superparticles. In this sum 
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contributions from reactions (|^) cancel, since they do not change the total number of 
superparticles. One thus has p| 



C^'^SUSY ou V^/ \ ( eq eq 



dt 
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-SHnsvsY - {(Tesv) (n^usY - ^susy) • (4) 



Here, H is the Hubble parameter, (...) denotes thermal averaging, v is the relative ve- 
locity between the two annihilating superparticles in their center-of-mass frame, and the 
superscript "eq" indicates the equilibrium density. In the second step we made use of the 
fact that, as argued above, all relevant heavier superparticles maintain relative equilib- 
rium to the neutralino LSP until long after the temperature Tp. This allowed us to sum 
all superparticle annihilation processes into an "effective" cross section; schematically 

(Teff cx + g^pBpaix'iP) + gpp {Bpf a(PP(*)). (5) 

where the Qij are multiplicity factors, and 

Bp = (mp/m,o)=^/2e-^'"^-™^?^/^ (6) 

is the temperature dependent relative Boltzmann factor between the P and Xi densities. 
The final LSP relic density VL^h?, where h = 0.65 ± 0.15 is the scaled Hubble constant, 
is then essentially inversely proportional to (cTcsv) at Tp ~ m^o/20. Co-annihilation can 
therefore reduce the LSP relic density by a large factor, if Sm = mp — m^o -C m^o and 

cr(xi-P) + cr{PP^*^) ^ cr(XiXi)- This is true in particular if x1 is a light, m^o < Mw, 
higgsino R H or SU(2)-gaugino W^. More recently it has been pointed out [ITT], O] that 



co-annihilation with light sleptons can reduce the relic density of a bino-like LSP by 
about one order of magnitude. 

In this letter we study co-annihilation of neutralinos with the lighter scalar top (stop) 
eigenstate ii. Compared to the other squarks, m^-^ is reduced by contributions of the 
large top quark Yukawa coupling to the relevant renormalisation group equations, as well 
as by mixing between SU(2) doublet and singlet stops. While we do not know of any 
model that predicts m^-^ ~ ""^x?' ^ close mass degeneracy is possible in many models, e.g. 
in the popular minimal Supergravity (mSUGRA) model 0. Moreover, scenarios with 
small ti-Xi mass splitting are of great concern for experimenters, since ti decays then 
release httle visible energy, making ii production very difficult to detect at both e"'"e~ [jl3[ 



and hadron |jTj] colliders. 

In contrast to the cases mentioned earlier, for P = ti it is not entirely obvious that 
reactions of the type (|^) will indeed be much faster than XiXi annihilation processes. 
In the absence of flavor mixing, one would have to chose X = W,Y = b or vice versa. 
However, for a temperature T < M^y, the H^-density is itself quite small, so reaction (||) 
would be much faster than XiXi annihilation only for m^o significantly above M^. On 
the other hand, most supersymmetric models predict some amount of flavor mixing in 
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the squark sector, even if it is absent at some high energy scale. As a result, for small Sm 
the dominant ti decay mode is usually its flavor changing 2-body decay into Xi + c . 



For ti masses of current experimental interest the dominant contribution to therefore 
comes from X = c,Y = nothing, i.e. (inverse) ti decay. If the effective ctix? coupling 
is suppressed by a small mixing angle e, the condition that (H) is much faster than XiXi 
annihilation reads 

^2g-WT.»^g-.o/T.^ (7) 

where the extra factor of a ~ 0.01 occurs since we are comparing 2^1 reactions with 
2^2 processes. For 6m ^ Tp m^o/20 we then only need e > e~^° ~ 5 ■ 10~^. 
In what follows we will assume that this is true, or that Xi is sufficiently heavy that 
Xi + ^ii + b is fast. 

Another property of the stop is that it has strong interactions. A leading order cal- 
culation of and cr(titi*^*'') will therefore not be very reliable. Unfortunately a full 
higher order calculation is highly nontrivial, since one would need to include finite temper- 
ature effects (e.g. in order to cancel Coulomb singularities in the non-relativistic limit). 
We expect these unknown higher order QCD corrections to be more important than the 
contributions of higher partial waves. In the calculation of the cross sections (j{xiti) 
and a{iiii'^) we therefore only include the leading, 5— wave contribution; however, the 
P— wave contributions to XiXi annihilation process are included. Our co-annihilation 
cross sections will thus only be accurate to a factor of 2 or so. Due to the exponential 
dependence of aes on 6m, see eqs.(|^-§|), the bounds on the ti-Xi mass splitting that will 
be inferred from upper or lower bounds on fi^/i^ should nevertheless be fairly accurate. 

The existence of unknown, but probably large, higher order corrections also means that 
we can ignore all ti annihilation reactions that involve more than the minimal required 
number of electroweak gauge couplings. However, we treat the top and bottom quarks 
Yukawa couplings on the same footing as the strong coupling (the latter Yukawa coupling 
will be large only for tan/3 ~ mt/mb). Altogether we therefore computed the cross sections 
for the following processes: 

xSti ^ tg,tH^,b H+; 
titi — > t t] 

Ul ^ gg, H''^ Hi H+ , bb,tt, (8) 

where Hf = h,H,A is one of the three neutral Higgs bosons of the minimal supersym- 
metric Standard Model (MSSM) |T6|. The cross sections for Xiti and tltl annihilation are 
identical to those in the first and second lines of eq.(^, respectively. We have performed 
two independent calculations of these cross sections. One calculation was based on trace 
techniques and the usual polarization sum for external gluons; here the non-relativistic 
limit (to extract the S"— wave contribution ||^) was only taken at the end. The second 
method uses helicity amplitudes in this case the non-relativistic limit can already be 
taken at the beginning of the calculation. [Note that the cross sections for titi — >■ g 
vanish in this limit.] Explicit expressions for these cross sections will be published else- 
where. [Note that in their analytical expressions, the authors of Refs.|n, T2\ do not keep 
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the mass of the relevant SM fermion /; these terms are irrelevant for their case f = t, 
whereas we have to keep a finite value for the top quark mass, 7^ 0. Ref.[jll| also did 
not include — fn mixing, which in our case is crucial for obtaining a light ti. In the 
relevant limit we agree with Ref.[jlT|. We disagree with Ref.[0 for the tg (or tj) final 
state, by a factor m^o/y/s.] 

Our calculation of the relic density closely follows ref.|Q]. In particular, we keep the 
temperature dependence (^) of (Jeff when computing the "annihilation integral" (essen- 
tially the integral of the Boltzmann equation for T > Tp). 

We use a variant of the minimal Supergravity model for our numerical analysis. 
In particular, we assume a common gaugino mass, a common sfermion mass mo, and a 
common trilinear soft breaking parameter Aq at the Grand Unification scale Mx = 2-10'^^ 
GeV. However, we allow the soft breaking masses of the two Higgs doublets to differ from 
mo- In practice, this means that we keep the higgsino mass parameter /i and the mass ttia 
of the CP-odd Higgs boson as free parameters at the weak scale. The final free parameter 
is the ratio tan/3 of vacuum expectation values of the two Higgs fields. 

For illustration, we take /i = — 2M2, where M2 ~ 2m-j^o is the SU(2) gaugino mass. 
This implies that the LSP is bino-like, which is the most natural choice for this type of 
model |T^. It is also conservative, since a higgsino-like LSP will have larger couplings to 
the (s)top, and hence even larger co-annihilation cross sections eq.^). We also chose a 
large sfermion mass, mo = 2M2. In the absence of co-annihilation this choice is usually 
incompatible with the upper bound on the LSP relic density, which we conservatively 
take as f^^^^ — 

In Fig. 1 we show contours of constant ^xh'^ in the (m^o, 6m) plane for a scenario with 
moderate tan/3 and a very heavy Higgs spectrum, mA = 5M2. This latter choice implies 
that the only Higgs boson relevant for the calculation of the LSP relic density is the light 
CP-even scalar h, with mass m/j < 130 GeV. This is a conservative scenario in the sense 
that it minimizes the number of final states contributing in eqs.(|^), and also leads to a 
small XiXi annihilation cross section. We see that scenarios with very large 6m are indeed 
excluded by the upper bound on Q^h"^. The peak in the contour fl^h'^ = 0.5 at m^o ~ mt 
is due to XiXi ~^ tt, which has a sizable 5— wave cross section if ti is not too heavy and 
m^o is not much above mt. The much smaller bumps at m^o ~ 130 GeV are due to hh 
final states becoming accessible. 

On the other hand, for very small values of 6m and m^o in the range indicated by 

naturalness arguments ( ^ 0.3 TeV, corresponding to a gluino mass mg ^ 2 TeV), we find 
that the LSP cannot contribute significantly to the solution of the Dark Matter puzzle, 
since its relic density is too small. In particular, one needs ^x^^ 0.025 for Xi to form 
galactic haloes. We see that even for the present very conservative choice of parameters 
one needs a ti-x? niass splitting of at least 9 to 19 GeV (6 to 10%) to satisfy this lower 
bound on VlyJi^. This mass splitting is large enough for standard ti search methods at 
e^e~ colliders ||13|, |19[ to have reasonably high efficiency. If we require that Vt^h"^ lies in 



the currently favored "best fit" range between about 0.1 and 0.2, 6m has to be between 
11 and 33 GeV. Unfortunately this is still not high enough for current ti search strategies 



at the Tevatron [p3 to be sensitive. 
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mg=-jU=2Mg, m^=5M3, tanj8=10 



100 
80 
^ 60 



> 
O 



1 1 1 






\ r)-u2_n 

\ WIl — U.iJ 




- 


\ : 

nh2=o.i \ 


, , 1 


nh2 = 0,025 '\ \ J 


100 


1000 




[GeV] 



^ 40 
20 




Figure 1: Contours of constant = 0.5 (solid), 0.1 (dashed) and 0.025 (dotted) in the 
{m^o,Sm) plane, where 6m = mi^ — m.^o. We took /i, mo and mA to be fixed multiples 
of M2 — 2m^o, as indicated, whereas tan/3 = 10 has been kept fixed. The parameter Ao 
varies between about 2.5mo and 3.2mo, with larger Aq values corresponding to smaller 
values of 6m. 



So far we have focused on LSP masses in the range favored by naturalness arguments. 
It is sometimes claimed |jl8| that the upper bound on fi^/i^ implies that LHC experiments 
must find superparticles if the MSSM is correct and Xi is bino-like. Unfortunately this 
is not true; for 6m ^ an LSP mass up to 4 TeV, corresponding to a gluino mass in 
excess of 20 TeV, cannot be excluded from this cosmological argument. [As noted above, 
our estimates for ti annihilation cross sections are not very reliable. However, even if 
we over-estimated them by a factor of 2, the bound on m^^o would only be reduced by a 

factor of ^/2, and would thus still allow sparticle masses far above the range to be covered 
by the LHC] 

In Fig. 2 we show analogous results for a light spectrum of Higgs bosons and large 
tan/3, where the bottom Yukawa coupling is sizable; the choice m^ = 0.35M2 ~ 0.7m^o 
ensures that all Higgs pair final states will be accessible for m^o > 100 GeV. However, we 
keep the previous (large) values for \fi\ and mo. Nevertheless we see that for natural values 
of m^o, requiring fi^/i^ > 0.025 now implies 6m > 20 GeV. Moreover, the LSP makes a 

good DM candidate, i.e. fixh'^ ~ 0.1, only for 6m ^ 40 GeV; this is sufficiently large to 
permit ti searches at the Tevatron [0, Finally, for 6m 0, cosmology now allows 
an LSP mass up to 6 TeV, corresponding to a gluino mass of about 30 TeV! Obviously 
the upper bound on 6m that follows from f^^^^ ^ 0.025 for natural values of m^o, as well 
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as the absolute upper bound on m^o that follows from ^^h"^ < 0.5, are even higher if we 
chose smaller values for tuq and/or 
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Figure 2: As in Fig. 1, except that we took a large value of tan/3 and a light Higgs boson 
spectrum. The regions below and to the left of the heavy dotted lines are excluded by 
Higgs boson searches at LEP. 



In conclusion, we have shown that scenarios with very small ti-Xi niass splitting would 
permit an LSP mass of several TeV without "overdosing" the Universe. This shows once 
again p that the upper bound on the LSP relic density does not guarantee that LHC 
experiments will detect superparticles, even if the MSSM is correct; of course, (third 
generation) superparticles with masses out of the reach of the LHC can hardly be argued 
to be "natural". On the other hand, for Xi and ti masses of present experimental interest, 
and indeed for the entire natural range of these masses, Xi cannot contribute significantly 
to the Dark Matter in the Universe unless the ti-Xi mass difference is large enough for 
conventional ti search strategies at e'^e~ colliders to be effective. This does not imply that 
collider searches for ti nearly degenerate with Xi |21| should not be continued; a positive 



signal would definitely exclude Xi DM candidate, which is not easy to accomplish with 
cosmological Dark Matter searches. However, since Dark Matter is known to exist, for 
natural values of m^o a very small ti-Xi mass splitting would require physics beyond the 
MSSM. 
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